Opportunities and limitations of transition voltage spectroscopy: a theoretical analysis 
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In molecular charge transport, transition voltage spectroscopy (TVS) holds the promise that 
molecular energy levels can be explored at bias voltages lower than required for resonant tunneling. 
We investigate the theoretical basis of this novel tool, using a generic model. In particular, we 
study the length dependence of the conducting frontier orbital and of the 'transition voltage' as a 
function of length. We show that this dependence is influenced by the amount of screening of the 
electrons in the molecule, which determines the voltage drop to be located at the contacts or across 
the entire molecule. We observe that the transition voltage depends significantly on the length, but 
that the ratio between the transition voltage and the conducting frontier orbital is approximately 
constant only in strongly screening (conjugated) molecules. Uncertainty about the screening within 
a molecule thus limits the predictive power of TVS. We furthermore argue that the relative length 
independence of the transition voltage for non-conjugated chains is due to strong localization of the 
frontier orbitals on the end groups ensuring binding of the rods to the metallic contacts. Finally, we 
investigate the characteristics of TVS in asymmetric molecular junctions. If a single level dominates 
the transport properties, TVS can provide a good estimate for both the level position and the degree 
of junction asymmetry. If more levels are involved the applicability of TVS becomes limited. 

PACS numbers: 



I. INTRODUCTION 

Molecular electronics aims at investigating and explor- 
ing the quantum properties of molecules in electronic 
devices [1-4]. An essential parameter when considering 
charge transport through molecules is the location of the 
molecular energy levels. In general, these levels are a few 
electron Volts (eV) away from the Fermi energy (Ep) of 
the electrodes. Hence, it should in principle be possible 
to align the Fermi energy of one of the electrodes with a 
molecular energy level by applying a bias voltage. How- 
ever, in practice a junction often breaks down before the 
molecular level is reached, as a result of the gigantic elec- 
tric field coming about (^ lO^V/m). 
Recently, Beebe et al. introduced transition voltage spec- 
troscopy (TVS) as an alternative method to characterize 
molecular energy levels in a device geometry [5, 6]. They 
determined current (/) vs. voltage {V) characteristics 
for a series of molecular devices, and replotted their data 
in a Fowler-Nordheim (FN) manner [7], i.e. they plot- 
ted ln(//l/^) versus 1/V. In such graphs a clear min- 
imum appears, at a voltage, Vmin, which is generally 
smaller than the voltage needed to reach the molecu- 
lar level (see Fig. 1 for a calculated example). Beebe 
et al. proposed that Vmin provides direct information 
about the energy distance between Ep and the nearest 
molecular level. This interpretation was based on a pic- 
ture of molecular junctions as tunnel barriers obeying 
the Simmons model for charge transport [8] . Experimen- 
tally, their proposition was supported by an extended se- 
ries of experiments [5, 6]. They demonstrated that Vmin 
does not vary with molecular length d for alkanethiols, 
which indeed have a length-independent HOMO-LUMO 
gap (The terms HOMO and LUMO refer to the highest 
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FIG. 1. (a) Typical current vs. the bias voltage curve for 
a model calculation (see main text). Once a molecular 
level is aligned to the Fermi energy of one of the electrodes, 
a step appears (at V = 1 V in this example where the 
molecular level energy is E = 0.5 , Ep = 0). (b) Fowler- 
Nordheim (FN) plot created from (a), showing ln(J/V^) 
vs. 1/V. A minimum appears at a voltage (Vmin), which is 
smaller than the voltage required to reach the molecular level. 



occupied molecular orbital and lowest unoccupied molec- 
ular orbital, respectively). Furthermore, they experimen- 
tally showed that for conjugated molecules Kiin scales 
linearly with \Ep — i^HOMO/LUMoL depending on which 
level, HOMO or LUMO, is closest to the Fermi energy. 

Recently, Huisman et al. found that the Simmons pic- 
ture is inconsistent with the experimental data of Beebe 
et al. [9]. They showed that within that model, Vmin 
is expected to decrease like 1 /d for a series of alkanethi- 
ols. (If the image potential is taken into account, this 
functional dependence on d changes for small d, see also 
Ref. [10]). Within a coherent molecular transport model, 
however V^in was found to be length independent for the 
same molecular series, provided d > 8 A [9]. Additionally, 
both Huisman and Araidai et al. pointed out that Vmin 
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does not occur at the transition between direct tunnel- 
ing and FN tunneling [9, 11]. Within a molecular level 
model, Knin rather appears when a certain amount of the 
tail of the broadened resonant level has come within the 
bias window. In two recent papers, the Thygesen group 
took the discussion a step further, by performing ab ini- 
tio calculations for a set of molecular junctions [12, 13]. 
In particular, they pointed out that junction asymmetry 
is an essential parameter for TVS to be interpreted cor- 
rectly [12]. 

TVS has clear potential in analyzing molecular charge 
transport experiments. However, the very issue if it 
can indeed be applied generally is still not settled. For 
this reason, we present an investigation of TVS using a 
generic theoretical model. Our analysis relies on a DFT- 
based many-body method described recently [14]. We 
focus on two essential questions that were brought up 
recently. The first one follows directly from the work of 
Huisman et al. [9]. An essential difference between the 
two models they compared is in the voltage profile as- 
sumed. In their coherent molecular level picture, the full 
voltage drops at the contacts, whereas in the Simmons- 
based model, the potential decreases linearly over the 
molecule itself. Hence, it is not clear whether the dis- 
tinction in length dependence of V^„in is due to the volt- 
age profile or due to the other clear differences between 
these models. Here, we study the influence of the exact 
voltage profile for a generic molecular model. The second 
issue we investigate involves the consequence of asymme- 
try for TVS (sec Rcf. [12]). Specifically, we introduce 
two separate T^nin values, for both positive and negative 
bias. Subsequently, we study how (and within which con- 
ditions) these quantities are related to both the position 
of the nearest molecular level(s) and the junction asym- 
metry itself. 

The organization of this paper is as follows. In Sec. II 
we introduce our model. As a generic system we use 
a Hubbard chain connected to two non-interacting con- 
ducting leads. To investigate the conductance through 
the molecule, we combine local spin density approxi- 
mation (LSDA) with many-body Green's functions. In 
Sec. HI A, the effect of a voltage drop over a molecule is 
studied, whereas the role of (a)symmetry is discussed in 
Sec. IIIB. 



II. MODEL AND METHOD 

The system we consider consists of a small region 
where possibly Coulomb interactions are present, weakly 
coupled to two non-interacting, semi-infinite leads (see 
Fig. 2). The interacting region contains one or several 
quantum dots in series. The Hamiltonian of the leads 
and the central region read respectively (H.c. denotes 
the Hermitian conjugate) 

ifloads = -ic 5Z HXl[4'),TCi-(-l.r,,<T +H.C.] (1) 
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FIG. 2. A short interacting Hubbard chain connected to two 
non-interacting leads. The gate voltage Vg can be applied to 
the interacting region. The hopping terms in the interacting 
part and in the contacts are t and tc respectively. The 
interacting region is coupled to the left and right contact by 
tL and tR. 
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where Nl is the length of the interacting chain. The 
parameters t and tc represent the hopping rate in the 
molecule and contacts respectively, and U describes the 
on-site Coulomb interaction. However, since we discuss 
the off-resonant regime and since in TVS the first step 
in the I-V characteristic is dominant, we can put U = 
without loss of generality. The creation and annihilation 
operators. c|, d|, Ci and di acting on site i satisfy the 
usual anti-commutation relations. In addition, the ex- 
ternal gate potential, Vg, can be applied to the central 
region which is included in the energy, i.e. e = e(V^). 
The index a =ti i describes the spin. 
The coupling Hamiltonian between the molecule and the 
contacts reads 

iJcoupling = E \-^'n4,ri,^dj^„ + H.c] (3) 

ri=L.R 
a 

where i denotes the leftmost (rightmost) site of the right 
(left) contact and j corresponds to the leftmost (right- 
most) site of the central region. 

Our method is based on a mapping of the Hamiltonian 
of the central region to a limited set of many-body eigen- 
states. All of the parameters in the many-body energy 
spectrum are obtained from ground state L(S)DA calcu- 
lations [15]. We then calculate the transport using many- 
body Green's function theory (see Ref. [14] for details). 
The retarded Green's function of a single level (without 
any Coulomb interaction) connected to the electrodes is 
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ular model of Ref. [9] , Vmin is found to be independent 
of the molecular length for d > 8A. Huisman et al. men- 
tioned that the differences in the functional dependence 
of Knin on thc length of the molecule could originate from 
the different voltage profiles in these two models. 
In this section, we investigate the length dependence of 
V^nin in more detail. Different factors influence the length 
dependence of transition voltage: screening, the hopping 
integral and the spatial structure of the conducting or- 
bitals. Here we systematically discuss thc influence of 
these factors. 



FIG. 3. Schematic representation of the energy levels of a 
short Hubbard chain connected to two non-interacting leads 
(with symmetrically applied bias voltage) considered in two 
different configurations, (a) without any voltage drop on the 
central region: the voltage drops only at contacts, (b) The 
potential decreases linearly with distance. 



described by 



(4) 



where Eg is the single electron energy [16]. The pres- 
ence of the contacts is taken into account by the retarded 
self-energy, E*". It is necessary to calculate the effective 
coupling of a level a to the contacts. We do this by 
projecting the central chain Hamiltonian onto two many- 
body states of N and 4- 1 or -t- 1 and N + 2 [14]. 
By extracting f^^, the calculation of the self-energies is 
straightforward. Once the Green's functions are known, 
the current can be calculated from a Landauer type of 
equation 
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where and G° are the lesser and advanced Green's 
functions respectively, Tj = i(E^ — E^ ) and f(uj,^j,j) is 
the Fermi distribution of lead j . 

In particular, for the transport through a single level or 
if the left and right line-width functions are proportional 
to each other, i.e. ~ AFjj, the current can be written 
as 



Li R 



(6) 



III. RESULTS AND DISCUSSION 
A. Length dependence of Vmin 

In the Simmons model, Vmin is found to be inversely 
proportional to thc molecular length while in the molec- 



1. Screening effect and hopping integral 

The capability of the electrons to screen out the field 
determines the way in which the voltage drops over the 
molecule. To investigate the difference between a voltage 
drop over the contacts and a linear drop over the entire 
molecule, we use our generic model and consider the two 
configurations shown in Fig. 3. Fig. 3(a) corresponds to 
the case of a molecular junction at nonzero bias with the 
drop entirely over the contacts. Fig. 3(b) corresponds to 
the case where the voltage decreases linearly with dis- 
tance between the leads. For model (b), a voltage depen- 
dent term is added to the diagonal elements of the Hamil- 
tonian. This term equals £:j?-T^/2 + (i-l)X^/(A^L-l) for 
site i, where is the number of central dots. In both 
cases (a) and (b), the hopping term in the central region 
is t. We note that in Rcf. [11] a comparable model is pre- 
sented. However, in that paper, two cases were studied 
in which not only the voltage drop was different, but also 
the hopping rate t and the energy values of the molec- 
ular levels were chosen differently, making it difficult to 
distinguish between the individual effects. Here we want 
to analyze thc behavior in a more general context by ad- 
ditionally calculating Vmin versus molecular length. 
First of all, for both models (a) and (b), we indeed see 
a minimum in the FN plot which is consistent with the 
results of Ref. [11], confirming that TVS should not be 
necessarily interpreted as a transition between rectangu- 
lar and trapezoidal barriers. The results of Vniin versus d 
for models (a) and (b) of Fig. 3 are shown in Fig. 4. This 
plot shows a very similar behavior for the two, rather dif- 
ferent, models (a) and (b). Both curves decay to become 
approximately constant only beyond a length of about 15 
dots [17]. 

It is useful to study, in addition to the d-dependence of 
Vnihi, the dependence on d of thc parameter 



X 



\Ef — Eg 



(7) 



suggested in Ref. [12]. This parameter gives an indi- 
cation of the relation between the TVS minimum and 
the nearest resonance, removing a possible length depen- 
dence of the latter from the problem (see below). 
We have plotted x versus d in Fig. 5. We see that there 
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FIG. 4. Vmin versus the length of the molecule (the number 
of the dots in the central region), t = 1, tc = 3, U = and 
Vg = —3. Squares for the model (a) of Fig. 3 and circles for 
the model (b) of Fig. 3. 



FIG. 5. X — \Ef — Esl/Vi-nin versus the length of the molecule 
(the number of the dots in the central region), t = 1, tc = 3, 
[7 = and Vg = —3. Squares for the model (a) of Fig. 3 and 
circles for the model (b) of Fig. 3. 



is a significant difference between length dependence of x 
for a voltage drop only at the contacts versus a drop over 
the entire molecule. It should be noted that the main 
differences between the curves in Fig. 4 (rather similar 
for model (a) and (b)) and Fig. 5 (which shows a striking 
difference between these two models) persist for different 
values of t and Vg. 

In fact the parameter x is a key quantity in our analysis 
as it gives more insight into the performance of TVS than 
Vmin- Moreover, if x is a well-defined number, the gap 
between the Fermi level and the nearest molecular level 
is easily calculated from Vmin- Returning to Fig. 5, we 
see that the parameter x shows the expected difference 
between voltage drop over the entire molecule versus a 
drop over the contacts only, whereas inspection of Vn-iin 
versus length, may not enable us to distinguish clearly 
between the two cases, as the HOMO or LUMO level it- 
self may vary significantly with length. 
Two factors influence the behavior of the quantities Vinin 
and X versus length: the extent to which the electrons 
screen out the applied potential, and the hopping inte- 
gral t. In Ref. [11] it was pointed out that the two are 
related: strong screening usually implies a large value of 
the hopping integral, as both result from a high mobility 
of the electrons. We discuss both effects in detail. 

The screening length in molecular systems varies 
strongly across different molecules due to the character- 
istics of different chemical bonds. For instance, screening 
is usually strong when there are many tt electrons in the 
molecule, in particular when they are arranged along con- 
jugated pathways [18, 19]. The amount of screening is re- 
lated to the HOMO-LUMO gap: A small HOMO-LUMO 
gap is indicative of strong screening [20]. Screening 
strongly influences the energy landscape through which 
the electrons move: in general if the electrons have the 
ability to screen out the voltage, the voltage drop will be 
localized near the contacts, while in the opposite limit 
the voltage drops over the entire molecule. As we have 
seen, the difference between strong and weak screening 
is most clearly observed in the dependence of x on the 



length, and very weakly in Vinin- 

On the other hand, the hopping integral t causes the 
HOMO and LUMO levels to vary with the chain length. 
In particular, for a chain of length d, the relation be- 
tween the frontier orbital and t can be seen in the energy 
spectrum given by 

E,, - Eo ^ -2t cos{ka) (8) 

where a is the inter-site distance (which is 1 in our case) 
and k = n-K / L where L = [d + l)a and Eq is the energy 
offset. Thus the maximum width of the energy spectrum 
is At, and the minimum energy for a chain of length d (if 
it were isolated) is given by 

Ei=Eo- 2t cos [n/{d + 1)] . (9) 

This level is relevant for LUMO transport; for 
HOMO transport the relevant level would be Eq + 
2t cos [7r/(c? -I- 1)]. We denote the energy level closest to 
the Fermi energy of the leads (the frontier orbital) by 
Es- As shown schematically in Fig. 6, when increasing 
the molecular length, the dominant transport level (i?s), 
moves closer to the Fermi energy (where in this exam- 
ple all of the levels are below the Fermi energy). Fig. 7 
shows the value Es of the frontier orbital with respect to 
the Fermi energy as a function of chain length. In fact, 
this variation of E^ with length, is the dominant factor 
in Fig. 4. 

To judge the effect of the hopping integral in a real 
molecule we consider known values for this parameter 
for the case of conjugated and non-conjugated molecules. 
The hopping integral t is reported to be 3.18 (eV) for a 
phenyl ring and 1.68 (eV) for alkyl [21]. Although these 
values show the expected trend, their difference does not 
seem dramatic enough to be responsible for the substan- 
tially different conductance of the two species. It seems 
therefore that the major difference between the conduc- 
tance properties of conjugated and non-conjugated sys- 
tems is not so much due to the difference in hopping inte- 
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FIG. 6. Schematic representation for variation of \Es — Ef\ 
versus the length of the molecule. By increasing the length 
of the molecule (the number of the dots in the central region) 
the closest level to Ef becomes closer to Ep and thus Vmin 
is decreased, t = 1, U = and Vn = —3. 
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FIG. 7. \Ef — Es\ versus the length of the molecule (the 
number of the dots in the central region) for the case that the 
voltage drops only at the contacts, t — 1, U = and Vg = —3. 
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FIG. 8. The hopping integral of an alkanethiol system is 
almost 1.5 eV [21]. Moreover DFT calculations show that 
the energy difference between frontier level of alkane located 
on sulphur and the rest of the alkane chain is about Ea = 1.5 
eV. Considering iJfronticr = 3.5 and Ea ~ 1.5, it is possible 
to calculate Eq from Eq. 8 for the alkane chain which is 
about 8 eV. Then the parameters Eq and t can be utilized 
to investigate the transport through such a chain using our 
generic model. Here have shown Vmin versus the length of 
the molecule for alkane chains with those parameters for (a) 
a Hubbard chain in which the first and last site are at a 
chemical potential of 3.5 eV and the rest of the chain is at 
the chemical potential of 8 eV and (b) a chain where all of 
the sites are at a chemical potential of 6.5 eV. t = 1.5, tc = 5 . 



gral. These two classes of molecules show however quite 
a different HOMO-LUMO gap, which indicates that the 
amount of screening is much larger in the conjugated 
case. 

We conclude that realistic values of the hopping integral 
lead to a Vmin which decays with length. This length de- 
pendence is expected to be noticeable in both conjugated 
and non-conjugated molecules. When looking at the pa- 
rameter Xi this length dependence disappears when the 
screening effects in the molecule are strong. In the next 
subsection we investigate why the length dependence of 
Vmin observed in the experiment is weak [6]. 



2. Spatial structure of frontier orbitals 

We have seen that V^nin varies with d tor d < 15, 
irrespective of the voltage drop. However, Beebe et al. 
concluded from an extensive analysis of experiments that 
Knin varies much less with d in alkanethiols [6]. This is 
also shown by H. Song et al. in Ref. [22] . We argue that 
the reason for this lies in the fact that the transport 
orbital is located mainly on the sulphur binding group 
[23-28]. Also, for the series of conjugated molecules used 



by Beebe et al. [5, 6], DFT calculations [29] suggest 
that the frontier orbitals or levels close in energy to the 
frontier orbitals are distributed over the entire molecule 
(in contrast to the alkanethiols where the transport 
orbitals are localized on the binding group). This may 
explain the length dependence of these molecules and 
why, for larger t, this length dependence should be 
stronger as HOMO/ LUMO levels vary stronger. 
To study the influence of the localization in alkanethiols, 
we include two sites which arc located at the ends of 
the molecular chain of our model. To be specific, we 
consider a Hubbard chain in which the first and last 
site are at a chemical potential of 3.5 eV and the rest 
of the chain is at the chemical potential of 8 cV and we 
compare this case with a chain where all of the sites have 
the same chemical potential of 6.5 eV. These numbers 
come from the DFT calculations for alkanethiols where 
the energy of the frontier orbital is about 3.5 eV and 
the energy difference between the frontier level of alkane 
located on sulphur and the rest of the alkane chain is 
about 1.5 eV [29]. Incorporating the energy level value 
with 2t = 3, Eq for the first case ( (a) in Fig. 8) would 
be about 8 eV and for the second case ( (b) in Fig. 8) is 
6.5. These numbers can also be seen from Eq. 8. The 



results are shown in Fig. 8. Knin for the first case is 
almost length-independent while for the second case, 
Knin varies significantly with d for d < 9. 
In a recent paper [30], it has been shown that a direct 
coupling (i.e. without a sulphur or other binding atom) 
between a carbon atom and gold is possible, with a high 
conductance when the carbon is sp3 hybridized. In that 
case, Vmin may vary more strongly with length. So far, 
TVS results for these alkanes have not been published. 
All in all, we conclude that for a homogeneous chain, 
the voltage drop has only a modest direct influence on 
the dependence of Knin on d. Comparing the variation 
of X = ^ function of d, we observe however 

a strong difference between a voltage drop over the 
contacts or across the molecule. For molecules with 
strong screening, the HOMO and LUMO are usually 
located close to the Fermi energy, and this increases 
the relative variation of Vmin with d. Finally, the weak 
length dependence of Vinin for non-conjugated molecules 
must be accounted for by the strong localization of the 
orbitals on the end group. 

Finally, from Fig. 5, we infer that for molecules with 
strong screening, x is a well-defined parameter which is 
constant within ±5% and allows us to infer the location 
of the HOMO (or LUMO) from the TVS minimum 
voltage. However, it is usually not well known where the 
voltage drops for a specific molecule, and this implies 
an uncertainty about the value of x- this value will lie 
somewhere between the upper curve and the lower curve 
of Fig. 5. Therefore this figure indicates the typical 
degree of uncertainty one faces in interpreting TVS 
results when the amount of screening is not known. 



B. (A) Symmetry 
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FIG. 9. X = \Ef — Ea\/Vmin versus rj for positive bias voltage 
and Es = 0.5. 
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We now turn to an investigation of the effect of asym- 
metric capacitive coupling on TVS. This aspect was pre- 
viously considered in Ref. [12]. Here we want to study 
this using our generic model. The asymmetry is de- 
scribed by a parameter ry, which we define differently 
from Ref. [12]. We change the chemical potential of 
the contacts by the parameter rj where jJ,L = Ep — rjV 
and /Xfl = Ep + (1 — ri)V [31]. In fact, the parameter 
77 specifies how the bias voltage is distributed over the 
left and right contact. In the symmetric case 77 = 1/2, 
while in the fully asymmetric case rj = 0. Using sym- 
metry 77 -f-^ 1 — 77, only 77 between and 0.5 needs to be 
considered. The ratio of \Ep — Es\ to Vmin, i.e. x versus 
77 for positive voltage is shown in Fig. 9. This ratio varies 
between 0.8 and 2, which is in agreement with the re- 
sult of Ref. [12]. Since x depends on symmetry, one can 
only find the HOMO/ LUMO level energy from Knin, if 
77 is known. In this curve, the molecular level is reached 
at V = \Ep — Es\/{1 — 77). Thus in order to have l^iin 
smaller than the voltage required to reach the level, we 
should have x > ^ ~ V which is the case for all of the 



FIG. 10. ln{I /V^) versus 1/V for positive and negative bias 
voltages in two cases of 77 = 0.1,0.5. Eg = 0.5, tL,R = 0.2 
and [7 = 0. 



points in Fig. 9. 

In Fig. 10, the FN plot is shown for a symmetric 
(77 = 0.5) and an asymmetric (77 = 0.1) junction. It can 
be seen that, for the asymmetric case, the transition volt- 
age differs between positive and negative voltages which 
is the first important feature that should be taken into 
account in the case of using TVS. 

In Fig. 11, X = \Ep — Es\/Vynin is shown as a func- 
tion of 7] for positive (xp) and for negative voltages 
(Xn) for two different cases : (i) Eg is above the Fermi 
level {Eg = 0.5, Ep ~ 0), i.e. the resonant level is 
the LUMO level, and (ii) Eg is below the Fermi level 
{Es -0.5, Ef=0), i.e. the HOMO level. Here we sup- 
pose that only a single level contributes to the transport. 
The absolute value of x is the same for the positive and 
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FIG. 11. X versus rj for Es = ±0.5 and Ef = 0. The left and 
right couphngs are II = tn = 0.2. Crosses: Es > 0, V > 0, 
diamonds: iSs > 0, V < 0, stars: Es < 0, V > 0, circles: 
Es<0,V <0. 



b 

3 
U 



^ IT 


= 0.5 




= 0.5 


ri 


= 0.2 




= 0.2 




Vbias 



FIG. 12. Current versus voltage for a two-levels model and 
Ei,2 — ±0.7, Ef ~ and Il.r = 0.15. Squares for positive 
voltage and rj — 0.5. Triangles for negative voltage and 
rj — 0.5. Crosses for positive voltage and 77 — 0.2. Circles for 
negative voltage and 77 — 0.2. 



the negative voltages in the symmetric junction as shown 
in Fig. 11. However it is different for the asymmetric case. 
Moreover, in the ease of the asymmetric jmiction, \x\ for 
> 0, F < is the same as for < 0, V > and also 
Ixl for E's > 0, 1/ > is the same as for Es < 0, V < 0. 

Interestingly, it is possible to estimate the asymme- 
try of the junction by looking at the ratio |Xp/Xn| = 
Vinin,n/Vmin,p whcrc the Subscript p, n refers to positive 
and negative voltages respectively. When Eg is above the 
Fermi level, the ratio of |xp/xn| or Vniin,n/Vniin,p is equal 
to (1 — r])/r]. Similarly in the case of Eg below the Fermi 
level, |Xp/Xn| is equal to 77/(1 — 77). Hence in principle, 
the ratio of 14iin,n/Vmin,p allows us to find the asymmetry 
degree and, hence remarkably, TVS can be used both for 
finding the information about the HOMO/LUMO level 
and about the asymmetry degree of the coupling by look- 
ing at positive and negative voltages. In Ref. [12], an 
asymmetrically coupled junction was studied for the pos- 
itive bias only and the authors argued that if the HOMO 
level can be found by other tools, then TVS can be used 
to estimate the asymmetry degree or vice versa which is 
different from our statement. 

Summarizing the results presented in this section, we 
emphasize that if one depicts the I-V characteristics by 
a FN plot, both the positive and negative bias voltages 
should be considered. From that, it is possible to say 
whether the capacitive coupling is symmetric or asym- 
metric. By having Vmin n/Vmin^p and knowing whether 
the resonant level is HOMO or LUMO, we showed that 
it is possible to estimate the degree of asymmetry of the 
molecular junction. However one should note that here 
also the voltage drop matters. Similar to the discussion 
in Sec. Ill A. the lack of information about the voltage 
drop can lead to an uncertainty in x of about ~ 20% to 
30%. 

In some cases, asymmetric coupling may lead to a sym- 
metric I-V. Consider for example a two-level system. The 



I-V curve for this system is shown in Fig. 12. The two 
closest energy levels to the Fermi energy {Ep = 0), are 
Ei^2 = ±0.7. In this case, even for asymmetric capacitive 
coupling (77 = 0.2), the I-V curve is symmetric and this 
is due to the fact that both energy levels are considered 
in the same distance far from the Fermi energy, with the 
same coupling. Hence the possibility to use TVS to es- 
timate the asymmetry degree of the molecular junction 
coupling only exists for the case where a single level is 
dominant in the transport. 

In experiments, usually one level dominates in the trans- 
port at low voltages and usually the HOMO and LUMO 
levels are not in the same distance from the Fermi en- 
ergy. In the case of a two level model, one could think 
of HOMO and LUMO energy levels below and above 
the Fermi energy with different distance from Ep. If 
\Ef - Eh\ > max{^,j^)\EF - El\ the transport 
is dominant through the LUMO level. A similar ar- 
gument leads to dominant HOMO when \Ef — El\ > 
max(i^, — Eh\- Both these conditions are 

the criteria to dominate the transport through one level 
which provide the possibility to use TVS. 



IV. CONCLUSIONS 

In conclusion, we have investigated the length de- 
pendence of transition voltage and the influences of 
(a)symmctric coupling of a molecular junction on TVS. 
For molecules with strong screening, the HOMO-LUMO 
gap is usually small compared to alkanethiols and this 
can explain the relative variation of V^nin with d in 
conjugated molecules. The weak length dependence of 
^min in alkanethiols can be elucidated by the strong 
localization of the orbitals on the end group. Unfor- 
tunately in the experiments the lack of information 
about the potential profile over the molecule limits the 
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usefulness of TVS as a tool for identifying the location 
of the conducting frontier orbitals. Furthermore we 
looked at the possibilities and shortcomings of TVS 
in the case of an asymmetric junction. We showed that 
in addition to get the information about the HOMO/ 
LUMO level position by TVS, it is also possible to 
estimate the degree of asymmetry by looking at the 
transition voltage value at the positive and negative 
voltages. It must be noted, however, that also this case 
is limited by lack of information on the voltage profile 
inside the molecule. Finally, we showed that in order to 
estimate the degree of asymmetry by TVS, one should 
take note of the number of the dominant levels in the 



transport. Summarizing, TVS may be a useful analysis 
technique, but it should be used with considerable care. 
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